Effect of bath temperature on the quantum decoherence 
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Abstract 

The dynamics of a qubit under the decoherence of a two level fluctuator (TLF) in addition to its 
coupling to a bosonic bath is investigated theoretically based on a unitary transformation. With 
the merit of our unitary transformation, non-Markov effect can be taken into account. It shows 
that quantum decoherence of the qubit can either be reduced or be increased with increasing bath 
temperature T or with increasing TLF-bath coupling. 
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Quantum computation reveals advantages over classical one for its highly efficient parallel 
calculation and therefore attracts a wide interest among scientistsjl, 2]. Today, the biggest 
obstacle is the quantum decoherence due to environmental noise, especially the intrinsic slow 
noise caused by the two level fluctuators (TLFs) . To overcome this difficulty, great effort 
has been made to design a well behaved qubit, and many methods have been proposed such 
as optimum working point J4], dynamical decoupling {5}], and etc. {(3, 7|. In this paper, we focus 
on the temperature dependent properties of a qubit under the influence of such noise. It is 
commonly believed that environmental disturbance can be reduced by decreasing the bath 
temperature and the system-bath coupling. However, an interesting paper shown that the 
decoherence can be reduced by increasing the thermal bath temperature for quantum systems 
coupled directly only to a few modes[8|. In that paper, the authors introduced a full quantum 
model in which the concerned two-state system (TSS-A) is only coupled to another TSS-B, 
with the latter being dissipated by a multi-mode thermal bath (See Fig. 1). Therefore, 
TSS-B and the multi-mode bath together can be seen as a non-Markov environment of 
TSS-A. In this paper, we treat the model more rigorously without making rotating-wave 
approximation (RWA) and Markov approximation, and find that the decoherence of TSS- 
A can either be reduced or increased with increasing bath temperature T depending the 
couplings conditions. When A-B coupling is smaller than B-bath coupling (condition 1), 
decoherence reduction with temperature happens as shown in Ref. 8|. When A-B coupling 
is larger than B-bath coupling (condition 2), only decoherence enhancing with temperature 
is expected. Further more, we find that, the decoherence of TSS-A can also be reduced with 
increasing B-bath coupling when condition 1 is satisfied. 

The model is given as (h = 1) 0, 0, Q, 11, Q, Q: H = H A + H AB + H B with 
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where TSSs are characterized by pseudospin-1/2 operators a A and as usual, and b\ 
are the annihilation and creation operators of the bath mode, go and gu are the coupling 



constants. Here we use the convention used in Ref. 



where only transverse coupling 



are included. The bath is fully defined by the spectral density J {of) = J^kdl^i 00 ~ u k)- 



2 



Suppose the TSSs are the double quantum dots (DQDs) manufactured with GaAs, we have 



the following spectral density [14 

J(uj 

As in Ref. 
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171 ]. we apply a unitary transformation to the Hamiltonian, H' = 
exp(S)H exp(— S), with the generator S = J2k ^^k(b[ — b^af. Note that this unitary 
transformation only affects Hb of the total Hamiltonian. After the transformation, Hb can 
be decomposed into three parts: H' B — H' Q + H[ + H' 2 , with 



H' 



m 
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^ J>(1 - + 6*) + ^i]A B l, 



-^A B (coshX -rj) + -^A B (sinhX - V X), 



(4) 
(5) 
(6) 



where X = ^£,k{b\ — b k ) and rj is the thermodynamic average of cosh AT, 

n = exp -J2Aticoth(f3uj k /2) , (7) 

which insures H' 2 contains only the terms of two-boson and multi-boson non-diagonal tran- 
sitions and its contribution to physical quantities is {gl) 2 and higher. Therefore, H' 2 will be 
omitted in the following discussion. 

The introduced parameter is an important point here. If = for all k, that is, 
without the transformation, the perturbation expansion would be similar to the standard 
weak-coupling expansion (Bloch-Redfield theory). Besides, if = 1 for all k, then our 
transformation is the usual polaronic transformation and the perturbation expansion is for 
the small parameter A which is equivalent to the NIBA 21]. Here, we choose £ k s as|22j|: 

OJk 
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UJ k + T]A B 



H[ is of rotating-wave form, 



H[ = J2Vk(b{a B + b k a*), 



(8) 



(9) 



with V k = g k r\A B j {yj k + r]A B ) and a± = (af =F«cr^)/2. It is easy to check that H[\g ) = 0, 
where \go) is the ground state of H' Q . Note that, the transformed Hamiltonian H' B = H' +H[ 



is of similar form of that of RWA, but A B and gfc/2 are renormalized by r]A B and Vk due 
to the contributions of anti-rotating terms. Here, r)A B gives a rough approximation of the 
renormalized qubit frequency and (1 — 77) A B is the corresponding Lamb shift of TSS-B in 
the sense of Quantum Optics 

At this point, we write out the master equation for this nonlinear system by treating 
H = Ha + H' B as unperturbed part, and Hab as perturbation. Suppose On -C A a, As, to 
the extent of second order of go, the reduced master equation of TSS-A is [19. I20]]: 



d PA (t) 
dt 



i[H A ,p A (t)}- [ dt'X(t,t'). (10) 
Jo 



with 



X(t,t') = Tt b [HAB,e- iHo ^[HAB,PA(t')®p B }e iHo ^ 
+ g 2 G 2 (t-t')e- iH ^ [p A (t%<T*]e iH ^*t 



where, we have assumed that the density matrix of A-B can be decoupled as Pa® Pb with 
p B = exp(— j3H' B ) /Tr B [exp(—f3H' B )] being the thermal equilibrium density matrix of TSS-B. 
Gi(t) and G 2 (t) are the correlation functions {(rf {t)af)'p and (o'fo'f(t))' /3 , in which (•••)« 
represents the average with thermodynamic probability ps- According to the fluctuation 
dissipation theory (FDT), G\(t) and G2{t) can be expressed by the Green's function of 
frequency ((erf ; erf ))' u . By solving the equation chain of motion, we can get the frequency 
domain of d(t) and G 2 (t) as G 1 {u) = G(u)/(1 + e~^) and G 2 ( u) = G{u)/{1 + eP") with 
(See Appendix) 

G{u) = lj# (n) 

where, | • ■ ■ | mean absolute value, R(u) and j(co) are the real and imaginary parts of 

T. t v t \2„ t + i)/(u,-m+-u, t ), 

J (u' + 7]A B y(uj - 10') 

7(0;) = 7r(77A B ) 2 J(cu)coth(^/2)/(cu + r7A i? ) 2 , (13) 

where p means the Cauchy principal value. 

The master equation Eq. fTTDT) which is a 2 x 2 matrix equation can be solved exactly by 
the Laplace transform method since the convolution theorem can be applied to the equation 



of each matrix element. Suppose the system is in the upper eigenstate of a z at the initial 

/io\ 

time t=0, p A (0) = • We can get the Laplace transform of population difference 

V 00 / 

(a?(t)) = Tr A (a?p A (t)) as 

P + 2F 



P 2 + 2PF + A 2 A 

where, F = g 2 Q f^du G{uj) / \P + iuj). After the inverse Laplace transformation and a 
parameter change of P = + + iu, we have the population difference P(t) = (oy(t)), 

Pit) - Re \2^ dU co 2 - 2uE F (co) -A\- 2iuT F (u) J ' (14) 

where, Re means the real part of a complex number, IV (a;) = T(u) + T(—u) and XV (w) = 
S(c<j) — E(— lo) are the imaginary and real parts of F(u) = g% J^duj' LU G ^_] + , with 



r(w) = vrj'H, e(w) = p / do/ i^l (15) 



where we have applied lj+ ^ 0+ = jp^ — Z7r5(a;), p means the Cauchy principal value, and 
J'{u)=glG{u)6{u). 

Before doing any numerical calculation, we would like to summarize the approximations 
we have made. Three approximations are made: The first one is the omission of H' 2 , which 
is a 4th order approximation to the B-bath coupling. H' 2 contains only multi-boson non- 
diagonal transition (like bkbk' and b\b\,) whose contribution is also small in low temperature. 
The second one is to approximate (b\ + bk)(b\, + by) ~ (271/. + 1)5^/ when calculating the 
correlator in the Appendix. The third one is the Born approximation for deriving the 
master equation (flQj) . Therefore, our treatment is applicable in low temperature for a^l 
and g <^A A ,A B . 

Now, we are in position of calculating P(t). We first use Eq. (J3j) , ( [T]) and (JH]) to get rj 
self consistently, then calculate P{t) numerically by using Eq. (l3l).(fnT)-(fl3l) and (fl~4 l) -( [TB I) . 
We report P(t) as a function of time in Fig. 2 and Fig. 3 for = A B = O.lui, g = 0.1A A 
and uJd = 0.05co7. In Fig. 2, a is larger than go/A A which is set to be a = 0.3, it shows that 
the decoherence is reduced by increasing the bath temperature T as predicted in Ref. 
However, in Fig. 3, where a = 0.01, the coherence is not meliorated but rather damaged 
with increasing T. In Fig. 4, we set a fixed detuning between A-B: 5 = A A — A B = 0.05a>z, 
the coherence is damaged with increasing T only because of the finite detuning 5 (compare 
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with Fig. 2). In Fig. 5, B-bath coupling constant a is varied with other parameter fixed, we 
can see the decoherence can also be reduced by increasing B-bath coupling. 

Now we are hoping for a bit more physical insight into the problem. For the sake of 
simplicity, we discuss it under the RWA to Hab, which is believed to be valid for the near 
resonance case (A^ pa Ag). If RWA is applied to the interacting Hamiltonian Hab, note that 
no RWA to B-bath interaction, according to the same procedure, the population difference 
P(t) would be, 

Pit) = Re{ / duo — — \ (16) 

K) \ 2ttz 7_ 00 uj-A A -E(iv)-iT(u)J 1 ; 

with S(co>) and r(a>) also defined by Eq. (TT5T) . One can find that, this expression has exact 
the same form of the population difference of SBM under RWA, only the spectral density 



J{uj) is replaced by J'{uj) = g 2 G(u)9(u) 



[l6, 



13, 



Therefore, from the point of view of 



the qubit, it feels like under the dissipation of a bosonic bath with an structured spectral 
density J' (id). If we assume the poles of G(u) to be ujb ± ij(u)B), with the renormalized 
frequency ojb being the solution of equation: uj — t]Ab — R(uj) = (Because of the property 
of the piezoelectric bath J(u), only one solution is expected and its value is always less than 
Ab). We can simplify J'(oj) as, 

J'(u) = 9hMIl , (17) 
[u — Ub] + 7 2 (^b) 

Note that this approximation would loss some non-Markov effects of the bath on the TLF. 

However, we will see that the main result of our numerical calculation can be explained, 

since the vital point is the non-Markov effect of the TLF+bath environment on the qubit. 

Approximately, the integral in Eq. (fT61) can be evaluated as P(t) = ^ a p e~ aplpt cos ou p t, 

p 

where, 7 P = T(u p ), a p is the weight of each pole and u p is the solution of equation 

uj — Aa — Ti(oS) = 0. (18) 

The above equation usually has two solutions u>±, and correspondingly, a p can be evaluated 
as a+ = Aa ~ w ~ and a_ 



UJ+— UJ— UJ+—LU- 



24|. In the limit of a pz 0, that is 7(0^3) — > 0, J' (us) 



gQ/n5(uj — ujb) , then u p can be solved as, u p = \ (A a + wb) ± a/ (A~a — ujb) 2 + 4g, 
Suppose the resonance condition = uj b ~ A B is satisfied, then oj p = A A ± go (also see 
;he inset of Fig. 3), which fully agree with the result of the simple Jaynes-Cummings model 



251 ] . In this case, the decoherent rate is 7a = 7 P = irJ'(u} p ). From Eq. ffTTj) and Eq. ([TBI . 



9o + 7 2 (Aa) 
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where, 7 (A A ) = 7o coth(/3A A /2) with 7o = n(r]A B ) 2 J(uj B )/(uj B + r]A B ) 2 « 7rJ(A A )/4. 
When T is large, we have 7 (A^) ^> (70? therefore, monotonic increasing function of the 
coherence with temperature is expected as in Ref. [8|. Nevertheless, when T is large, other 
decoherence mechanism and multi-phonon process may take place, and our treatment is 
no longer suitable. Therefore, we are interested in whether temperature can indeed reduce 
decoherence at low temperature. When T is small, go and 7 (A a) would compete with 
each other in the denominator of Eq. ffl9l) . Which one would dominate depends on the 
coupling condition, e.g., 7 (A^) would dominate when g <C cx pz Aa and g dominate when 
o. pz Aa <C go- Therefore, for small T, ja can be expressed as 

T f 70 co4L A /2) ' 9o « a pz A A « A A « A B ^ 
I 7o coth(/3 A A /2) , A A < #0 < A A « A B 

That is, decoherence reduction with increasing T and a pz happens only when g <C a pz A A <^ 
Aa ~ Ab, which is the region Markov approximation applies|26[]. When a pz A A <^ go non- 
Markov effect becomes important, and different temperature dependent behavior comes out. 

In conclusion, without making RWA and Markov approximation, the dynamics of a qubit 
coupled with nonlinear bath was investigated through a perturbation method based on a 
unitary transformation. We find that: 1. The decoherence of TSS-A can either be reduced 
or increased with increasing bath temperature T. The couplings condition controls the 
temperature dependent bebavior. Decoherence reduction with T happens only when the 
A-B coupling is smaller than B-bath coupling. 2. Finite detuning between TSS-A and 
TSS-B would destroy this decoherence reduction effect. It would benefit by setting A a a 
little smaller than A B , since the effective level spacing of TSS-B uo B is always smaller than 
original A B because of the dressing of phonons. 3. Keeping g a p2 A^ A^ ~ A B but 
increasing a pz also reduces the decoherence of TSS-A. 
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APPENDIX A: THE CALCULATION OF CORRELATION FUNCTIONS 



For the sake of simplicity, we can apply a unitary transformation to H' B with the generator 

S=ina B /4 + u k b\b k , the Hamiltonian H' B becomes, 
k 

H'k = + £ Vkblb k + M b l°- + b ^+)> ( A1 ) 

k k 

where, the definition of a± becomes a± = (erf ± ia B )/2. Correspondingly, Gi(t) and G^i) 
become (erf (t)crf and (of erf (f))£. 
The equation of motion is given by 

B))l = ([A, B] + )» + (([A, H' B ]\ B))l (A2) 

where, ((A\ B))" represents the Fourier transformation of the Green's function 
—i9(t){[A,B] + )'L Then, we can get the equation chain of {{^\ a x))'Li (foyl a x))i 
(M&l + b k )\a*)) and ((a z (b\ - b k )\o»)). Since, [of i^of] = (2b\b k + and 
[a B b k ,b\a B ] = -{2b\b k + l)a B , we have, 

+ h), ib\a B + 6^)] = (2n k + l)ia*, (A3) 
fcf(&I " 6*), (b[a* + b k a B )] = (2n k + l)af , (A4) 

where, we have substituted b\b k by its thermodynamic average value n k and omitted all the 
b\b\ and b k b k terms. This approximation makes the equation chain self-closed. By solving 
these equations, we get 

uj-r]A- V, v * {2nk+1 > u + V A - V. v " (2 " fc+1) 
According to FDT. Gi(t) and G2{t) can be expressed by the corresponding retarded 
Green's function ((erf | &x))'ui+io+-> so we § e ^ 

GiM = G( W )/(l + e>), (A6) 
G 2 (u,) = G(u)/(l + e*>), (A7) 

with = -l/7rlm((erf | a B ))' u;+i0+ which leads to Eq. ffTT]) . 
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FIGURES CAPTIONS 

Fig. 1: diagrammatic sketch of a qubit coupled with a two level fluctuator 

Fig. 2: P(t) as a function of time in the case of g <C a pz A A , where the decoherence is 
reduced with T. Inset is the Fourier analysis of the main plot, where the half width of each 
plot characterizes the rate of decoherence. Since the integration of each plot P(ou) must be 
one as a result of the initial condition, the higher the peak is, the narrower it should be, 
therefore, the half width which indicate decoherence is also smaller. We can see the inset is 
consist with the main plot. 

Fig. 3: P(t) as a function of time in the case of g ^> a pz A A , where the decoherence is 
enhanced with T. Inset: From the Fourier analysis of P(t) one sees that two frequencies 
are dominating the dynamics and the peaks locate at ± g which agree with the result 
of Jaynes-Cummings model. 

Fig. 4: Fourier analysis of P(t) in the case of g <C a pz A A with a finite detuning. The 
time domain of P(t) is similar to the plots in Fig. 2, we only give the frequency domain 
here for clear and simplicity. The decoherence is enhanced with T just because of the finite 
detuning compared to Fig. 2. 

Fig. 5: Fourier analysis of P(t) with different B-bath coupling a pz in the case of 
9o ^ a P z^A, the decoherence is reduced with increasing a pz . 
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